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1 Introduction 

Several recent papers have studied the possibihty of baryogenesis in the Standard Model and 
its minimal extensions. For example, references [1] - [6] impose constraints on the parameters 
of the Standard Model, minimal supersymmetric model (MSUSY), models with additional 
bosons or singlet Majoron model. In particular, authors of ref. [1] estimate the upper limit 
on the lighter CP-even Higgs in MSUSY to be equal to that in the Standard Model which 
they calculate to be 55 GeV. However, they have restricted their analysis to the case of the 
top quark lighter than about 115 GeV. The present paper includes the corrections due to the 
heavy top and its supersymmetric partners. 

The obvious question is: are these corrections important? Contributions of the heavy top 
and stop to the effective potential have been calculated in many papers ( [7] - [13] ) and turn 
out to be a dominant part of the one loop effective potential . These results have been used 
in ref. [14] in order to impose limits on the parameters of MSUSY from existing data at LEP 
and in ref. [15] to analyze Higgs signals in the future hadron supercolliders. 

In this paper we consider the one loop effective potential including contributions from W 
and Z gauge bosons, top and stop quarks and two CP-even Higgs bosons of the MSUSY in 
the region of parameter space where they play a non-negligible role. 

Next we numerically calculate the zero and nonzero temperature VEV's and impose the 
requirement that the latter be large enough immediately after the phase transition to suffi- 
ciently suppress baryon number violating processes. This constraint gives upper bounds on 
the lighter CP-even Higgs mass and on the soft SUSY breaking scale /i. 

In section 2, we briefly review how baryogenesis imposes constraints on parameters of the 
Standard Model and in section 3 we give a short description of the Higgs sector of MSUSY. 
Our calculation is presented in section 4 and the results are analyzed in section 5. 

2 Constraints from Electroweak Baryogenesis 

Much attention has been devoted recently to the possibility of creating the baryon asymmetry 
during the electroweak phase transition. For this to be true one has to satisfy certain basic 
requirements. First, the amount of CP violation inherent in the model has to be large enough 
to account for the observed baryon to photon ratio of about 10~^. Standard Model does not 
satisfy this requirement but it is possible to extend it to have enough CP violation. Secondly, 
the phase transition has to be first order to provide for departure from thermal equilibrium. 

Once these requirements are met it is possible to construct a mechanism which shows 
how the baryon asymmetry was created during the electroweak phase transition, ( see for 
example ref. [3] ) instead of at much higher energies. This lower energy scale will allow for 
predictions of baryogenesis to be experimentally verifiable. 

Regardless of the details of the particular mechanism of creating the baryon asymmetry, it 
is important to make sure that baryon number violating processes after the phase transition 
do not erase any previously created symmetry. To be more specific: after the phase transition 
at temperature T^.., the Higgs field acquires a vacuum expectation value v{Tc). The rate of 
baryon number violation by thermal fiuctuations is proportional to the Boltzmann factor 
exp{—Msph/Tc), where Msph is the mass of the sphaleron field configuration or equivalently the 



height of the barrier separating the gauge field configurations with different baryon numbers. 
^sph was calculated at zero T in ref. [16]: 

Msph{T,)=AiTB{\)— (1) 

with: -B(^) G [1.52, 2.70], for \ E [0, oo), where the A is Higgs self couphng. 

These baryon number violating processes have to proceed at a rate much smaller than the 
expansion rate of the universe after the transition, therefore the exponent of the Boltzmann 
factor has to be sufficiently large. Shaposhnikov [17] showed that: 

'P^^ "" > 45 (2) 

This equation was derived under the assumption [22] that at T = Tc (1) can be approxi- 
mated by the same expression with v, A, g^ evaluated at T = Tc . 

How does this constrain parameters of the model? For instance, in the Standard Model, 
the one loop effective potential V'^{(f)) of the Higgs field at a temperature T much higher 
than the masses of the particles is given by: 

V^{<P) = 7(T^ - T^)<P' - ET<p' + ^0^ (3) 

where the coefficients 7, E, At are positive and are determined by the parameters of the 
Standard Model. In particular, A^ is the temperature dependent effective quartic self coupling 
of the Higgs field. 

The phase transition occurs near the point where the temperature dependent effective 
mass of the Higgs field vanishes: 

7(T2-T|)^0 (4) 

Then, the VEV is given by: 

v{n) ^ ^ (5) 

If we want to have v{Tc) sufficiently large to satisfy (1) and (2), we must have At suffi- 
ciently small. This in turn imposes the upper limit on the Higgs self coupling A and therefore 
on the Higgs mass. 

The similar thing will happen in MSUSY as our results will show. 

3 Higgs Sector of MSUSY 

Higgs sector of MSUSY (see for instance ref. [18]) contains two complex Higgs doublets with 
the following SU{3) x SU{2) x U{1) quantum numbers: 

H,= (^\e (1, 2, -1/2) , ^2 = ( ^i ) G (1, 2, +1/2) (6) 



From these eight real fields spontaneous symmetry breaking decouples three unphysical 
Goldstone bosons and one is left with five physical Higgs bosons, namely: two CP-even 
scalars, one CP-odd scalar and a pair of charged scalars. The tree level Higgs potential: 

V = m\\Hi\^ + ml\H2\^ -ml{HiH2 + h.c.) + ^{H+aHi + H+aH2f 
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+ ^-^{\H,\'-\H2\'f (7) 

o 

can be restricted to the real components of the neutral Higgs fields, 0i = Re Hi, 
h = ReH^: 
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One can always choose such a field basis that m\, fi, V2 are real and positive. Constants 

mi, 1712 and m^ have to satisfy certain conditions. Requiring that the potential be bounded 

from below gives: 

.2 I ^2 

> ml (9) 



"^1 + "^2 ^ 2 



2 
and the spontaneous symmetry breaking condition is: 

m\ > m\m\ (10) 

Here Vi and V2 are proportional to vacuum expectation values of 0i and 02 : 

<^.>.^,<^.>.^ (U) 

and tan[3 is defined to be their ratio: tanjd = f2/f 1. Here: 



/vf + vi = 246GeV (12) 
to reproduce the measured values of gauge boson masses: 

2 9l / 2 , 2\ 2 9l ' 92 / 2 , 2\ /io\ 

my, = ^(^^1 + V2) , m^ = — - — [v^ + V2) (13) 

Fields 01 and 02 couple to down and up type quarks respectively. For example, after the 
spontaneous symmetry breaking top quark gets the mass: 

m? = /i? < 02 >2 (14) 

However, if the supersymmetry is softly broken, its scalar superpartner "stop" will have 
the different mass: 

mj = h] < 02 >^ +Ai^ (15) 

where we consider only the case of common soft supersymmetry-breaking mass /i for ti and 

Ir and vanishing off-diagonal elements of the 2x2 stop mass matrix. (The following analysis 
can be easily generalized to include all of these terms.) 



Finally CP-even physical eigenstates li and H, with masses rrih < mn are obtained by 
diagonalizing the mass matrix: 



M.- 



1 Id^V, 



tree 



A C 
C B 



and their masses are given by: 



ml^ = -{A + BT^J{A-BY + AC^) 



(16) 



(17) 



All this was at the tree level, but, as was already mentioned, in order to obtain limits on 
the Higgs mass it is paramount to include one -loop corrections. In the effective potential 
approach used in ref. [7] and [8] , masses of the Higgs bosons are approximated with the 
eigenvalues of the matrix of second derivatives of the one-loop effective potential evaluated 
at its minimum. 

The one loop effective potential at zero temperature is given by the expression: 



V" = Vtree{Q) 



1 

647r2 



Str I M^{(p) 



In 



Q' 



li 



Here A^^(0)is the field dependent squared mass matrix, Q is the renormalization scale 
and the supertrace is given by: 



Str f (m') = J2{-iy'^gJ (m 



(19) 



The sum runs over all the physical particles i of spin Jj, field dependent mass eigenvalue 
TTii and multiplicity gi that couple to fields 0i and 02- In our case these are W and Z bosons, 
top and stop quarks and h and H bosons with multiplicities: 
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6 , ^2 = 3, gt=g^ = 12 , gy^ = gH = l 



(20) 



We have neglected the contributions due to other quark-squark flavors. This is justifiable 
insofar as their masses are small. As was pointed out in ref. [7], also the bottom-sbottom 
contributions can be non negligible for very large values of tan(3. This case will not be 
relevant for us. 



4 Nonzero Temperature Effective Potential 

When the temperature is nonzero, effective potential gets a contribution [19]: 

mi(0) 



A\/. 



T 



2^ 



J29il 



i-fi 



T 



(21) 



where g^ are multiplicities as before, 'mi{(f)) are field dependent masses and /_(/_|_) which are 
to be used for bosons (fermions) are given by: 



I^{y) = ^ r x^ In (l ± e-v^^^+^') dx 



f22) 



This contribution to the effective potential describes the interactions of the Higgs bosons with 
the thermal bath surrounding them. Expressions (22) are rather difficult to operate with, 
especially when the masses depend on fields in a complicated way. Fortunately, as shown in 
ref. [2], one can always use either high temperature (small y) or low temperature expansion 
(high y), so that the mistake in determining AV^ is never bigger than 10 percent. 
High temperature expansions of (4.2) are given by: 

^-(^) = -45 + T2^ - 6^ -32^H^J' "^'^^"^^' ^ UJ' 

lncfe^5.41,lnc/^2.64, (23) 

Whereas the low temperature expansion is: 

,(,).-;|,3/V»(l + |) (24) 

In this calculation we will always use one of these expansions or linear interpolation 
between them. By substituting field dependent values of m^, ruz, rrit, mi, rrih and mn from 
equations (13) - (17) into the expressions for the effective potential (18) and (21) one obtains 
the full zero and nonzero temperature one-loop effective potentials. The critical temperature 
in this system is close to the point where the temperature dependent effective mass matrix 
has a zero eigenvalue. 

What can we get out of this? If we take rrit and fi to be our input parameters we have: 
3 SUSY parameters mi, m2 and m,^; 2 zero-temperature VEV's < 0i > and < 02 >, 2 
nonzero-temperature VEV's < 0i >t and < 02 >t and the critical temperature - altogether 
eight unknowns. 

How many conditions do we have? First the fixed magnitude of the zero T VEV (11) and 
(12), then 2 zero-T minima: 

^ = f^ = (25) 



next 2 nonzero-T minima: 



the critical temperature condition: 



dV^ dV^ 



(26) 



and, finally, from (1) and (2), the condition that nonzero-T VEV be sufficiently large: 

v{n) = ^vl{T^)+vl{T,) > vUTe) = jBlx/g^) ^^^^ 

If we take the equality in (28) (which corresponds to the upper limit on the m,h ), we have 
seven equations in eight unknowns, therefore we can impose one relation between them. This 
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In the MSUSY case: 



therefore: 



was done numerically in the form: rrih = rrihitanjS) for different values of parameters rrit and 
/i. Here, rrih is the upper limit on the mass of the lighter CP-even Higgs field. 

This program was realized in MATHEMATICA and maxima of curves mh{tanj3) are given 
in figures 1,2 and 3 for rrit = 115 GeV, 150 GeV and 200 GeV respectively. 

The argument \^^^ j g^ in the function B[\^^^ /g"^) was determined in ref. [5] for the general 
case of a two doublet model: 

\^ff = Ai cos^ pT + As sin^ f^r + 2h cos^ ^t sin^ /3t (29) 

Here (3t is the nonzero-T "mixing angle of VEV's": 

tan/^T = -j7^ 30 

A, = A, = £?±9i,ft = -£l±i (31) 

A=" = ^ii^cos2(2/?T) (32) 

There are several causes of uncertainty in this calculation. First, we have assumed that 
the phase transition happens at point T^ where (27) is satisfied. This is not true. It was 
already noticed in ref. [1] that since the phase transition happens earlier, when the vacuum 
expectation value is smaller than at T^, the actual bound is stronger than the one we take. 
In other words, if we were able to calculate the phase transition temperature exactly the 
upper limit on the Higgs mass would be lower than the one we impose. Unfortunately, at the 
phase transition the one loop effective potential is not accurate at the origin due to infrared 
divergences and therefore we can just estimate the critical temperature. 

Secondly, when using eq. (28) for the MSUSY we used the fact derived in [5] that the 
upper bound on the sphaleron mass in MSUSY is the sphaleron mass of the Standard Model. 
Therefore, our bound is again weaker than the actual one but still it will turn out to be very 
strong. However, the more precise calculation would require calculating the sphaleron mass 
in MSUSY at the critical temperature. 

Finally, as is usual in the study of phase transitions in early universe, we are using effective 
potential which is a static quantity for a system which not only evolves but evolves out of 
equilibrium. 

One should keep all of these caveats in mind when interpreting the results which are given 
in the next section. 



5 Results and Discussion 

Figures 1, 2 and 3 show upper limits on the Higgs mass for mt=115, 150 and 200 GeV 
respectively. Points are obtained as maxima of curves mh{tanf3) for different values of /i. For 
easier visibility they have been connected by straight lines. One can draw two conclusions 
from these results. 



First, for values of /i = 150GeV (which is the asymptotic lower experimental mass limit 
at 90% c.l. for a gluino mass lower than 400 GeV - see ref. [20] ), we get the upper limit on 
the Higgs mass to be 51 GeV, 54-5 GeV and 63 GeV for 3 different values of the top mass. 
One should compare this with experimental lower limit of 41 GeV^ (see [14] and [21] ). 

Second, for the considered region of tan(3 > 1 , there was always a maximal value of fi 
above which v{Tc) was never big enough to satisfy requirement (28). This gives the upper 
limit on the soft-SUSY breaking scale of about 750 GeV, 250 GeV and 170 GeV for three 
top masses considered. One should compare this with the previously mentioned asymptotic 
lower mass limit of 150 GeV. 

As a conclusion one can establish the following "no loose theorem" from these results: 
either the top quark is heavy (fig 3.) in which case Higgs can be as heavy as 63 GeV but 
SUSY breaking scale is very close to its experimental lower limit or top is lighter than 150 
GeV (fig 1. and 2.) but then the Higgs mass is lighter than 55 Gev and thus close to its 
experimental lower limit. 

As with all other calculations in supersymmetric models this one has the trouble that 
there are simply too many unknown parameters. We have considered that region in param- 
eter space which has been searched by experiments [14] ,[20], [21] ( i.e. for a common soft 
supersymmetry-breaking mass for ti and tu and vanishing off-diagonal mass elements for 
the stop mass matrix ). The limits that we obtain complement the experimental results and 
severely limit the parameter space of MSUSY. 
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